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Abstract 

The theory of fully developed turbulence is usually considered in an idealized 
homogeneous and isotropic state. Real turbulent flows exhibit the effects of 
anisotropic forcing. The analysis of correlation functions and structure func- 
tions in isotropic and anisotropic situations is facilitated and made rational 
when performed in terms of the irreducible representations of the relevant 
symmetry group which is the group of all rotations SO(3). In this paper 
we firstly consider the needed general theory and explain why we expect dif- 
ferent (universal) scaling exponents in the different sectors of the symmetry 
group. We exemplify the theory context of isotropic turbulence (for third 
order tensorial structure functions) and in weakly anisotropic turbulence (for 
the second order structure function) . The utility of the resulting expressions 
for the analysis of experimental data is demonstrated in the context of high 
Reynolds number measurements of turbulence in the atmosphere. 



I. INTRODUCTION 

Experiments in fluid turbulence are usually limited to the measurement of the velocity 
field at one single spatial point as a function of time. This situation has begun to improve 
recently, but still much of the analysis of the statistical properties of Navier- Stokes turbulence 
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|1[ is influenced by this tradition: the Taylor hypothesis is used to justify the identification 
of velocity signals at different times with differences of longitudinal velocity components 
across a spatial length scale R. Most of the available statistical information is therefore 
about properties of longitudinal two-point differences of the Eulerian velocity field and their 
moments, termed structure functions: 

S n (R) = (\[u(r + R) - u(r)\ ■ , (1) 

where (■) denotes ensemble averaging. In isotropic homogeneous turbulence, these struc- 
ture functions are expected to behave as a power-law in R, S n (R) ~ R^ n , with apparently 
universal scaling exponents ( n . 

Recent research f|]|| has pointed out the advantages of considering not only the longitu- 
dinal structure functions, but tensorial multi-point correlations of velocity field differences 

w(r,r',t)=u(r',t)-u(r,t), (2) 

given by 

^n "' S ( r h r[, h; r 2 , r' 2 , t 2 ; • • • ; r n , r' n , t n ) 
= (w a ( ri , r[, t!K(r 2 , r' 2 , t 2 ) . . . w 5 (r n , r' n , t n )) , (3) 

where all the coordinates are distinct. Even when the coordinates fuse to yield time- 
independent structure functions depending on one separation only, these are tensorial quan- 
tities [|]] denoted as 

S a(3 -(R) = (K(r + R) - u a (r)} [u p {r + R) - u p {r)} . . .) . (4) 

Needless to say, the tensorial information is lost in the usual measurements leading to ([[]). 
One of the main points of the present paper is that keeping the tensorial information can 
help significantly in disentangling different scaling contributions to the statistical objects, 
contributions that are hard to distinguish when quantities like (|I|) are considered. Espe- 
cially when anisotropy implies different tensorial components with possible different scaling 
exponents characterizing them, careful control of the various contributions is called for. 
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To understand why irreducible representations of the symmetry group may have an im- 
portant role in determining the form of correlation functions, we need to discuss the equations 
of motion which they satisfy. We shall show that the isotropy of the Navier-Stokes equation 
and the incompressibility condition implies the isotropy of the hierarchical equations which 
the correlation functions satisfy. We will use this symmetry to show that every component 
of the general solution with a definite behavior under rotations (i.e., components of a defi- 
nite irreducible representation of the SO (3) group) has to satisfy these equations by itself - 
independently of components with different behavior under rotations. This "foliation" of the 
hierarchical equations may possibly lead to different scaling exponents for each component 
of the general solution which belong to a different SO (3) irreducible representation. 

In Sect. 2 we describe the general mathematical framework of theory by discussing the 
structure of tensorial fields from the point of view of 5*0(3) irreducible representations. 
We then show in Sect. 3 that the hierarchy equations are indeed isotropic and as a result 
foliate into different sectors of the SO (3) irreducible representations. In the next sections 
we demonstrate the utility of the theory. In Sect. 4 we revisit Kolmogorov's four fifth's 
law emphasizing the rule of the SO (3) irreducible representations in its derivation. Then, 
in Sect. 5, we present some experimental evidences for the importance of an anisotropic 
exponent in the second order structure function, in atmospheric measurements. Sect. 6 offers 
conclusions and some comments about the road ahead. 

II. TENSORIAL CORRELATION FUNCTIONS AND SO (3) IRREDUCIBLE 
REPRESENTATIONS: GENERAL THEORY 

The physical objects that we deal with are the moments of the velocity field at differ- 
ent space-time locations. In this section we suggest a way of decomposing these objects 
into components with a definite behavior under rotations. We will show later that compo- 
nents with different behavior under rotation are subject to different dynamical equations, 
and therefore, possibly, scale differently. Essentially, we are about to describe the tenso- 
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rial generalization of the well-known procedure of decomposing a scalar function \P(r) into 
components of different irreducible representations using the spherical harmonics: 

tt(r) = $>, m (r)y, m (r) . (5) 

l,m 

A. Formal definition 

Consider a typical moment of the velocity field, Eq.(H). F^ 1 - an (ri,r' l 
is a function of 2n spatial variables and n temporal variables. Physically, it is a tensor field: 
if F n is measured in two frames / and / which are connected by the spatial transformation 
(say, a rotation) 

x° = k a px p (6) 
then, the measured quantities in each frame will be connected by the relation: 

Fn ( r l) r 'l' ' ' ' i r ni r ' n i (7) 

= A" 1 ^ • ■ ■A a ^ n F^(r l ,v' l ,t l] . . .;r n ,v' n ,t n ) 

= A^ ■ ■■A a ^ n F^(A- l r 1 ,A-%,t l ; . . . ; A _1 r n , A _1 f^, t n ) . 

We see that as we move from one frame to another, the functional form of the tensor field 
changes. We want to classify the different tensor fields according to the change in their 
functional form as we make that move. We can omit the time variables from our discussion 
since under rotation they merely serve as parameters. 

Consider coordinate transformations which are pure rotations. For such transformations 
we may simplify the discussion by separating the dependence on the amplitude of from 
the dependence on the directionality of r,: 

T Q — ( ri ,...,r p ) 
= T Q - a "(r 1 ,...,r p ;r 1 ,...,r p ) 
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For pure rotations we may treat the amplitudes ri, . . . , r p as parameters: the transformations 
properties of T ai ~ Mn under rotation are determined only by the dependence of T ai - an on 
the unit vectors ii,...,r p . Accordingly it seems worthwhile to discuss tensor fields which 
are functions of the unit vectors only. Notice that in the scalar case we follow the same 
procedure: by restricting our attention to scalar functions that depend only on the unit 
vector f , we construct the spherical harmonics. These functions are defined such that each 
one of them has unique transformation properties under rotations. We then represent the 
most general scalar function as a linear combination of the spherical harmonics with r- 
dependent coefficients, see Eq. (||). 

The classification of the tensor fields T Ql ' Qn (r 1 , . . . , r p ) according to their functional 
change under rotations follows immediately from group representation theory |5]||. But in 
order to demonstrate that, we must first make some formal definitions. We define 5™ to 
be the space of all smooth tensor fields of rank n which depend on p unit vectors. This is 
obviously a linear space of infinite dimension. With each rotation A G 50(3), we may now 
associate a linear transformation 0\ on that space via the relation (|7|): 

[0 A Tr>-* (f 1? ...,r p ) 

Using this definition, it is easy to see that the set of linear operators 0\ furnishes a repre- 
sentation of the rotation group £0(3) since they satisfy the relations: 

General group theoretical considerations imply that it is possible to decompose S™ into 
subspaces which are invariant to the action of all the group operators 0\. Moreover, we can 
choose these subspaces to be irreducible in the sense that they will not contain any invariant 
subspace themselves (excluding themselves and the trivial subspace of the zero tensor field). 
For the 50(3) group each of these subspaces is conventionally characterized by an integer 
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j = 0, 1, 2, . . . and is of dimension 2j + 1 [0,0. It should be noted that unlike the scalar case, 
in the general space <S™ , there might be more than one subspace for each given value of j. 
We therefore use the index q to distinguish subspaces with the same j. For each irreducible 
subspace (q,j) we can now choose a basis with 2j + 1 components labeled by the index m: 



J-'qjm \ L 1) • • • > L p) i 



m= j. 



In each subspace (q,j), the group operators 0\ furnish a 2 j + 1 dimensional irreducible 
representation of £0(3). Using the basis ^•m"' a "(ri, • • • , r p ), we can represent each operator 
0\ as a (2j + 1) x (2j + 1) matrix £>^ m (A) via the relation: 



[OaS]^-°» (r 1; ...,r p ) 



rrr =— j 



It is conventional to choose the basis B gJm such that the matrices D^, m (4>), that correspond 
to rotations of radians around the 3 axis, will be diagonal, and given by: 



The «S" space possesses a very natural inner-product: 



(T, U) = Jdki ...d% 



■T a - a - (xx... x p )^ Ql/3l . . . g anPn U 01 " 9n ( Xl . . . x p )* 



where g a/ 3 is the 3-dimensional Euclidean metric tensor: 



V 



/ 



By definition, the rotation matrices A a /3 preserve this metric, and therefore it is easy to see 
that for each A G SO (3) we get: 
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(0 A T,0 A U> = (T,U) 

so that, 0\ are unitary operators. If we now choose the basis B g j m to be orthonormal with 
respect to the inner-product defined above, then the matrices D^) m (A) will be unitary. 

Finally, consider isotropic tensor fields. An isotropic tensor field is a tensor field which 
preserves its functional form under any arbitrary rotation of the coordinate system. In 
other words, it is a tensor field which is invariant to the action of all operators (9 A . The one 
dimensional subspace spanned by this tensor-field is therefore invariant under all operators 
0\ and therefore it must be a j = subspace. 

Once the basis B gjm has been selected, we may expand any arbitrary tensor field 
F ai - a "(ri, . . . , r p ) in this basis. As mentioned above, for each fixed set of amplitudes 
ri, ... , r p , we can regard the tensor field F ai -- an (ri, . . . , r p ) as a tensor field which depends 
only on the unit vectors ?!,..., f p , and hence belongs to <S™ . We can therefore expand it in 
terms of the basis tensor fields B gjm with coefficients that depend on 7*1, . . . , r p : 

F Q1 - Q "( ri ,...,r p ) 

= £ ^i.",*^!.---,?,) • (8) 

q,j,m 

The goal of the following sections is to demonstrate the utility of such expansions for the 
study of the scaling properties of the correlation functions. 



B. Construction of the basis tensors 

a. The Clebsch-Gordan machinery . 

A straightforward (although somewhat impractical) way to construct the basis tensors 
'Bqjm is to use the well-known Clebsch-Gordan machinery. In this approach we consider 
the S p space as a direct product space of n 3-dimensional Euclidean vector spaces with p 
infinite dimensional spaces of single-variable continuous functions on the unit sphere. In 
other words, we notice that S™ is given by: 

r 1 1 n r nip 



and therefore every tensor T ai ~ cin (r 1 . . . r p ) can be represented as a linear combination of 
tensors of the form: 



v"* are constant Euclidean vectors and ^(r?) are continuous functions over the unit sphere. 
The 3-dimensional Euclidean vector space, Sq, contains exactly one irreducible representa- 
tion of SO (3) - the j — 1 representation - while «S°, the space of continuous functions over 
the unit sphere, contains every irreducible representation exactly once. The statement that 
Sp is a direct product space may now be written in a group representation notation as: 

p times 



n times 



<S™ =1 <g> 1 ® . . . <g> 1 <g> (0 © 1 © 2 . . .) <g> . . . (0 © 1 © 2 . . .) 



We can now choose an appropriate basis for each space in the product: 



• For the 3-dimensional Euclidean space we may choose: 



ei 



1 

72 



e 



V2 



V 



V 1 / 



e_i 



1 

72 



/ 1 ^ 



—2 



For the space of continuous functions over the unit sphere we may choose the well- 
known spherical harmonic functions. 



Once these bases have been chosen, we can construct a direct-product basis for 5^: 



jjjai...a„ 

ii...i n (hfii)...(l P iip) 



(?!,..., r p ) 



The unitary matrix that connects the Ei 1 ...i n (i 1 ^ 1 )...(i p ^ p ) basis to the ~Q q j m basis can be calcu- 
lated using the appropriate Clebsch-Gordan coefficients. The calculation is straightforward 
but very long and tedious. However, the above analysis enables us to count and classify the 
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different irreducible representations of a given j. By using the standard rules of angular- 
momentum addition: 

s ® I = \s - l\ e • • • © (s + l) 

we can count the number of irreducible representations of a given j. For example, consider 
the space Sf of second-rank tensors with one variable over the unit sphere. Using the 
angular-momentum addition rules we see: 

Si = 1 <g> 1 (0 © 1 © 2 © 3 © . . .) (9) 
= (0©1©2)©(0©1©2©3©...) 
= (3 x 0) © (7 x 1) © (9 x 2) © (9 x 3) © .. . 

We see that there are exactly three j = representations, seven j = 1 representations and 9 
representations for each j > 1. It can be further argued that the symmetry properties of the 
basis tensors with respect to their indices come from the 1©1 = 0©1©2 part of the direct 
product (|9]). Therefore, out of the 9 irreducible representation of j > 1, 5 will be symmetric 
and traceless, 3 will be anti-symmetric and 1 will be tracefull and diagonal. Similarly, the 
parity of the resulting tensors (with respect to the single variable) can be calculated. 

Once we know how many irreducible representations of each j are found in S™, we can 
construct them "by hand", in some other, more practical method which will be demonstrated 
in the following subsection. 

b. Alternative derivation of the B gjm . 

The method we wish to propose in this subsection is based on the simple idea that con- 
tractions with r a , 5 a @, e a ^ 7 and differentiation with respect to r a are all isotropic operations. 
Isotropic in the sense that the resulting expression will have the same transformation prop- 
erties under rotation as the expression we started with. The proof of the last statement 
follows directly from the transformation properties of r a , <5 Q/3 , e af31 . 

The construction of all B gjm that belongs to <S" now becomes a rather trivial task. We 
begin by defining a scalar tensor field with a definite j, m. An obvious choice will be the 
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well-known spherical harmonics Y} m (r), but a better one will be: 

$ im (r) = r%„(r). 

The reason that we prefer $ jm (r) to Y} m (r), is that $j m (r) is polynomial in r (while Y^ m (r) 
is polynomial in f) and therefore it is easier to differentiate it with respect to r. Once we 
have defined $j m (r), we can construct the ~B q j m by "adding indices" to ^^(r) using the 
isotropic operations mentioned above. For example, we may now construct: 

• r-^$ jm (r), 

• r-i +2 5 a Pd T d^ jm {r), 

• r~j~ 1 x a Qj m (r), etc... 

Notice that we should always multiply the resulting expression with an appropriate power 
of r, in order to make it r-independent, and thus a function of r only. 

The crucial role of the Clebsch-Gordan analysis is to tell us how many representations 
from each type we should come up with. First, it tells us the highest power of f in each 
representation, and then it can also give us the symmetry properties of B g j m with respect 
to their indices. For example, consider the irreducible representations of Sf - second rank 
tensors which depend on one unit vector f . The Clebsch-Gordan analysis shows us that this 
space contains the following irreducible representations: 

sl = sl ® si <8> si 

= 1<8>1<8>(0©1©2©3©...) 

= (0©1©2)®(0©1©2©3©...) 

= (3 x 0) © (7 x 1) © (9 x 2) © (9 x 3) © .. . 

That is, for each j > 1 we're going to have 9 irreducible representations. The indices 
symmetry of the tensor comes from the Sl (E) Sl = 1©1 = 0©1©2 part of the direct 
product. This is a direct product of two Euclidean spaces, so its a second rank constant 
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tensor. We can mark the representation number in this space with the letter s, and the 
representation number of the <S° = 0©1©2©3©... space with the letter /. This way each 
representation in Sf of a given j will have two additional numbers (s, I), which actually serve 
as the index q that distinguishes irreducible representations of the same j. The s index will 
determine the indices symmetry of the tensor, while the / index will determine the highest 
power of f in the tensor. If we now recall that in the space of constant second-rank tensors, 
Sq ©Sq = 0©1©2, the s = 0, 2 representations are symmetric while the s = 1 representation 
is anti- symmetric, we can easily construct the B ?jm : 



(s,l) = (0,j) B ljm (r) 


= r-^$ jm (r), 






(s,l) = (1, j - 1) B 2jrn (r) 


= r~ j+1 e 


^^(r), 




(s,0 = (l,j) B 3jm (r) 










(s,l) = (1, j + 1) B 4jm (r) 


= r~ j ~ 




^%$ jm (r), 




(s,l) = (2,j-2) B 5jm (r) 


ee r-^'+ 2 d Q d^ m (r) 


i 


(10) 


(s,l) = (2, j - 1) B 6jm (r) 


= r -j+i 






(s,/) = (2,j) B 7jm (r) 


ee r _J ' 


r <*Q/3 _|_ r /3^a 






(s,l) = (2,j + l) B 8jm (r) 


ee r- j - 


1 






^,j'm( r )j 


(s,l) = (2,j + 2) B 9jm (r) 


^ — J — 




a r^ ]m {v). 






It should be stressed that these B 9jm are not exactly the same one we would have gotten from 



the Clebsch-Gordan machinery. For example, they are not orthogonal among themselves 
for the same values of j, to (although, they are orthogonal for different values of j or to). 
Nevertheless, they are linearly independent and thus span a given (j, to) sector in the 
space. 

III. THE ISOTROPY OF THE HIERARCHY OF EQUATIONS AND ITS 

CONSEQUENCES 

In this section we derive the equations of motion for the statistical averages of the velocity 
and pressure fields differences. We start from the Navier-Stokes equations, and show that 
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its isotropy implies the isotropy of the equations for the statistical objects. Finally, we 
demonstrate the foliation of these equations to different sectors of j, m. 

Consider a Navier- Stokes incompressible turbulence in a bounded domain Q. The equa- 
tions of motion describing the flow are: 

d t u a + u^u" = -d a p + ud 2 u a , 
d a u a = . 

As is well known, the relevant dynamical time scales are revealed only when the effect of 
sweeping is removed. In our work we use the Belinicher-L'vov transformation |7j in which 
the flow is observed from the point of view of one specific fluid particle which is located at 

r at time t . Let p(r ,t \t) be the particle's translation at time t: 

t 

p(r , t \t) = J dsu[r + p(r , t \s), s] . (11) 
to 

We then redefine the velocity and pressure fields to be those seen from an inertial frame 
whose origin sits at the current particle's position: 

v(r , t \r, t) = u[r + p(r , t \t),t] , 
7r(r , £ |r, t) = p[r + p(r , t \t),t] . 

Next, we define the differences of these fields: 

W Q (r , t \r, r', t) = v a {r , t \r, t) - ^(r , t \v\ t) , 
n(r ,t |r, r',t) = 7r(r , t \r, t) - 7r(r , t \r', t) . 

A straightforward calculation shows that the dynamical equations for W are: 

d t W a (r, r', t) = - (d a + d' a ) n(r , t |r, r', t) (12) 
+v(d 2 + d' 2 ) W a (r ,t \r,r',t) 
-^W^r , t |r, r , t) W a (r , t |r, r', t) 
-d'^(r , t \r', r , t) W Q (r , t \r, r', t) , 
d a W a (r ,t \r,r',t) = d' a W a (r ,t \r,r' ,t) =0 . 
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By inspection, to merely serves as a parameter, and therefore we will not denote it explicitly 
in the following discussion. Also, in order to make the equations easier to understand, let 
us introduce some shorthand notation for the variables (r] t ,r' k ,tk)' 

Xfe = (r fc , r fe) t k ) , 
Xk = (rk, r' k , t k ) , 
X fe = (r fc , f fe ) . 

Using (|T2D, we can now derive the dynamical equations for the statistical moments of 
W,ri: Let (•) denote a suitable ensemble averaging. We define two types of statistical 
moments: 

•^'n 1 '" Qn ( r o|Xi, . . . , X n ) 

= (W ai (r |X 1 )...W a »(r |X n )> , 

7^" 2 '" Qn (r |Xi, . . . , X n ) 
= (nCrolXOW^CrolX^.-.V^CrolXn)) . 

Equation (|12D implies: 

a tl ^-^(ro|X 1; ...,X n ) (13) 

= " ( d (r\) + 9$)) HT" an (r \Xu • • • , X n ) 
-efpjvR--** (r |X,X 1; ...,X n ) 
-d { < ] m\- an (rolX'.Xi,...^) 
+1/ (d 2 {ri) + 9 ( 2 r , } ) ^-""(rolX!, . . . ,X n ) , 
XEE(r ,r',t) ; X'ee (r, r , t) , (14) 

C l) ^ Ql - Q ' l (ro|X 1 ,...,X n ) = 0, (15) 
«9ff ) ^ Ql - Q ' l (r |X 1 ,...,X n ) = 0. 

Equations (|T3"D, ( TT5] ) are linear and homogeneous. Therefore their solutions form a linear 
space. The most general solution to these equations is given by a linear combination of 
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a suitable basis of the solutions space. To construct a specific solution, we must use the 
boundary conditions in order to set the linear weights of the basis solutions. We shall 
now show that the isotropy of these equations implies that our basis of solutions can be 
constructed such that every solution will have a definite behavior under rotations (that is, 
definite j and m - see Sect. 2). But before we do that, note that in many aspects the 
situation described here is similar to the well-known problem of Laplace equation in a closed 
domain Q: 

V 2 ^ = , 
^\an = o . 

The Laplace equation is linear, homogeneous and isotropic. Therefore its solutions form a 
linear space. One possible basis for this space is: 

*i, m (r) = r l Y lm (r) , 

in which the solutions have a definite behavior under rotations (belong to an irreducible rep- 
resentation of SO (3) ). The general solution of the problem is given as a linear combination 
of the # l>m (r): 

tf(r) = 5>J.m*j,m(r) . 

l,m 

For a specific problem, we use the value of \l/(r) on the boundary (i.e., we use cr(r)) in order 
to set the values of aim- 



To see that the same thing happens with the hierarchy equations ( |13|, |15l) , we consider 
an arbitrary solution {^F n ,7i n \ n = 2,3,...} of these equations. According to Sect. 2 we 
may write the tensor fields J- n , 7i n in terms of a basis B gjm : 

^-^(rolXx,...,^) (16) 

qjm 



R (n) ai...a»/£ 
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7^" 2 - an (r |Xi, . . . , X n ) 

E HqjL{ r 0; Xl, ■■■iXn 



Xn qjm l r 0) -Vb ■ ■ ■ ^nj • 

Now all we have to show is that the pieces of J- n , 7i n with definite j, m solve the hierarchy 
equations by themselves - independently of pieces with different j, m. The proof of the last 
statement is straightforward though somewhat tedious. We therefore only sketch it in general 
lines. The isotropy of the hierarchy equations implies that pieces of ^F n ,1-C n with definite 
j, m, maintain their transformations properties under rotation after the linear and isotropic 
operations of the equation have been performed. For example, if jF" 1,Qn (r |Xi, . . . , X n ) 
belongs to the irreducible representation (j,m), then so will the tensor fields: 

$n)f*i-a* t Bfaj**-"*, etc... 

although, they may belong to different S™ spaces (i.e., have one less or one more indices). 
Therefore, if we choose the bases {b^} to be orthonormal, plug the expansion (|16D into 



the hierarchy equations equations ([13|, |TJ), and take the inner product with B^, we will 



-,(n) rj(n) _ 



obtain new equations for the scalar functions F Qj ^ H t 

iV 0, ^1, • • • , 

r( n ) U Y Y M3( n_1 ) u( n ) 



d tl F£J n (r ,X 1 ,...,X n ) (17) 



E ([Pin) + ^W)J H q ijm( r 0, Xl, ■ ■ ■ , X n )B q ,j m , 'Bgj,,,, 
q' 

~ E (d»^ F tfjm ( r °> X,X U ..., X n )B q ^\ B q n jr] 

q' 

~ E (d^F^U X', X U ..., X n )B^\ Bg, 
+^ E ((^(n) + ^K)) -^i'iL( r o ; • • • , X n )B£ jm , B q f m ^ , 



q' 



E (C^i't K X 1? . . . , X n )Bg m , B%?) = • 
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Note that in the above equations, (•) denote the inner-product in the S™ spaces. Also, the 
sums over q',j', m! from ( |T6| ) was reduced to a sum over q' only - due to the isotropy. We 
thus see explicitly from (|H],|I8|) the decoupling of the equations for different j, m. 

A. Rescaling Symmetry and Anomalous Exponents 

The hierarchical equations simplify somewhat in the limit of infinite Reynolds number 
Re— > oo. This limit is equivalent to v —* 0, in which the last term in Eqs.(|j]) can be 
neglected with impunity. It was pointed out before || that this is the main advantage of 
using "fully unfused" correlation functions in which all the coordinates are distinct: there 
is nothing to compensate for the vanishing of the viscosity in the v — ► limit. Once the 
viscous term is discarded, the rest of the equations . exhibit invariance to rescaling under 
the following rescaling group: 

1 1 /vi j, b l r /\ ^1 , J- qj m T /\ J- qj m , 

rr(n) x (n+l)h+Z 3 (h) rj-(n) ,. q s 

±± qjm /y ± - L qjm i V u I 

as can be verified by direct substitution. In ( |19D A and h are arbitrary scalars, and Zj{h) is 
an arbitrary n- independent scalar function. We endow it with an index j since we expect, 
and see below, that Zj(h) will differ in different j-sectors, but not in different m-sectors. 

As a consequence of the rescaling symmetry we can seek solutions that do not mix values 
of h. We define F^^ and H^n h as ^ ne quantities that solve the equations of motion on 
an /i-slice, which are the same as equations ( |I?D without the viscous term. The important 
property of the solution on an /i-slice is that it is a homogeneous function of all its arguments 
in the sense that 

F q jm,h(^r , Ari, Xr[, A 1_/l ti, . . . Ar n , Xr' n , A 1 "^,) = 

X nh+Zj{h) F qjm , h (r , n, r[, h,... r n , r' n , t n ) . (20) 

It should be stressed that the quantity F^ itself is not homogeneous in its arguments. It 
has been discovered in || and stressed in that time-correlation functions in turbulence 
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do not satisfy dynamic scaling in the sense of Eq. fl20|) . Indeed, the solution of Eq.(|17D is a 
sum over contributions on /i-slices, 



i#i(r , X 1 ...,X n )= dfi(h)F^ h (r , X u . . . , X n ) , (21) 

" ^min 

with n{h) some unknown measure that needs to be obtained from boundary conditions. 
Eq. ( pTD can be endowed with further meaning by rescaling coordinates and times according 
to 

Pj = r j/ R n, P'j = r'j/Rn, Tj = tj/tn^h , (22) 

where R n and ti? ni h are defined as the typical scale of separation of the set of coordinates 
and the typical times scale on that scale on an h -slice: 

^^Eh-r;f, (23) 
R (L\ h 



'^uKli) ' < 24 > 

Here U is the typical velocity on the outer scale of turbulence L. Defining now 

E j = (pj^Pj^ T j) > (25) 
Eq.(f2~T|) can be written, using the rescaling property on an h slice as 

Fqj m {X\, • • • , X n ) = 

/imax / B \ nh + Z i( h ) „. . 



U n dp(h) F^fr, ...,E n ) (26) 

This form is known as the "multi-fractal" form P-pH . The scaling exponents characterizing 
a q n ^ m are obtained from a saddle-point calculation in the limit R/L — > as mm.h{nh+Zj(h)}. 

It was explained in [|3] that Zj{h) is obtained from a solvability condition of the hierarchy 
of equations fll7|) . In particular the numerical value of the function Zj(h), and consequently 
of the scaling exponents which are determined by the saddle point integral, depend on the 
coefficients in the equations (|17D . We found that the scalar functions associated with the 



different j -irreducible representations, F^^Xi, . . . , X n ), satisfy equations with different 
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coefficients, depending on inner products of the basis functions B g j m . Accordingly we ex- 
pect the scalar function Zj(h) to change from sector to sector. If the functions Fjj^ are 
characterized by anomalous exponents, they may be different for different j. On the other 
hand, for the same j the equations mix different m (and q) components, and unless there is 
an additional symmetry to SO (3), we do not expect different contributions with the same 
j to exhibit different exponents. In the next section we will demonstrate explicitly in the 
context of 3rd order correlation functions how the existence of an additional symmetry, in 
that case parity, brings about a foliation of a j sector into two sub- sectors which exhibit two 
different scaling exponents. 

IV. EXAMPLE: KOLMOGOROV'S "FOUR-FIFTH LAW" AND THE FOLIATION 

TO DIFFERENT J'S 

One of the best known results in the statistical theory of turbulence is Kolmogorov's 
"four-fifth law" which was discovered in 1941 |I2"|] . This law pertains to the third order 
moment of longitudinal velocity differences 5ui(r, R, t) = [u(r + R, t) — u(r, t)] ■ H/R where 
u(r, t) is the Eulerian velocity field of the turbulent fluid. The fourth-fifth law states that 
in homogeneous, isotropic and stationary turbulence, in the limit of vanishing kinematic 
viscosity v — > 



where e is the mean energy flux per unit time and mass e = v (| Vq/U/3| 2 ). The only assumption 
needed to derive this law is that the dissipation is finite in the limit v — > 0. 

In this section we revisit this law by finding the full tensorial form of the j = component 
of the 3rd order correlation function. We do not need to assume isotropy of the turbulence 
at any stage of the development; the isotropy of the equations of motion suffices to decouple 
the j = contribution from all the rest, and in this case we have enough equations to 
determine the j = component of the tensor completely. We will also show that the j = 




(27) 
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component has two subcomponents with different scaling exponents. These subcomponents 
have different parity, and therefore are further decoupled in the equations of motion. The 
usual fourth-fifth law pertains to the components that have odd parity. One can derive an 



additional exact relation that pertains to the even parity components pB |. 

Defining the velocity v(r, t) as v(r, t) = u(r, t) — (u) we consider the simultaneous 3rd 
order tensor correlation function which depends on two space points: 

J a '^(R) = (v a {r + R, t)v?{r, t)v\r, t)) . (28) 

We show [[HJ that in the limit v — > 0, under the same assumption leading to the fourth- fifth 
law, this correlation function reads 

j*.07(r) = _ -L(Rrs af} + RPs<*-y - ^R a S^) (29) 

- ^-(e ap5 W + e^R^Rg . (30) 
The quantity h is the mean dissipation of helicity per unit mass and time 

h = u((d a u f3 )(d a [Vxu] f3 )) , (31) 



The new result fl30|) (derived firstly in [T3|] and [|Tj]) can be also displayed in a form that 
depends on h alone by introducing the longitudinal and transverse parts of u: the longitu- 
dinal part is uz = R(u-R)/i? 2 and the transverse part is u t = u — u^. In addition we have 
5uz(r,R, t) = 5ui(r, R, t)Ti/R. In terms of these quantities we can propose a "two fifteenth 
law" that pertains to the j = component of the following correlation function: 

([<fu,(r, R, t)] ■ [u*(R + r, t) x u t (r, *)]) = ±hR 2 . (32) 

We note that this result holds also when we replace u by v everywhere. 

To derive the result ([!(]) we start from the correlation function J a '^ 7 (R) which is sym- 
metric with respect to exchange of the indices (3 and 7 as is clear from the definition. Using 
the symmetry the most general form of the j = component of this tensor can be written 
by observation (Sect. 2b): 
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J a '^(R) = a^R^PR 1 + 5 ai RP + 8^R a ] (33) 

+a 2 ( J R)[<P /3 i? 7 + S^B? - 25^R a ] 

+az(R)[5 aP k 1 + 8 a ~<R? + <^ 7 ir - hR a Rf 3 R? / R 2 } 

+a 4 (R) [e a ^W + e" 7 ^^ (34) 

This form is precisely of the type J2 q a qjmBqjm for the isotropic sector j — m — 0. 

Not all the coefficients are independent for incompressible flows. Requiring 
dJ a '^(R)/dR a = leads to two relations amone; the coefficients: 

(m + ^^m - ^ + (35) 

(^ + |)[5«i(B)-4a 2 (fl)] = . 
As we have two conditions relating the three coefficients ai, a 2 , only one of them is 



independent. Kolmogorov's derivation [|T2| related the rate of energy dissipation to the 
value of the remaining unknown. Here the coefficient remains undetermined by the 
incompressibility constraint. It belongs to a component of odd parity; since the equations 
of motion and the incompressibility constraint are invariant under parity transformation it 
decouples altogether and needs to be determined separately. 

Kolmogorov's derivation can be paraphrased in a simple manner. Begin with the second 
order structure function which is related to the energy of R -scale motions 

S 2 (R) = (|u(R + r,t) -u(r,t)| 2 ) . (36) 

Computing the rate of change of this (time-independent) function from the Navier-Stokes 
equations we find 

= = -V 2 (R) - 2e + uV 2 S 2 (R) , (37) 

where T> 2 (R) stems from the nonlinear term (u- V)u and as a result it consists of a correlation 
function including a velocity derivative. The conservation of energy allows the derivative to 
be taken outside the correlation function: 
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V 2 (R) = ^(/(r,tK(r + R ) t)[/(r,t) (38) 
-V(r + R,t)]) . 

In terms of the function of Eq. (|28"D we can write 

V 2 {R) = ^ [j"^(R, t) - J a ^(-R, t)} . (39) 

Note that Eq. (^) is written in terms of v rather than u, but using the incompressibility 
constraint we can easily prove that Eq. (|38| ) can also be identically written in terms of v 
rather than u. We proceed using Eq. (^) in Eq. (^), and find 

V 2 (R) = 2-^R^ ai {R) + 25i(i2)] . (40) 



For R in the inertial interval, and for v — > 0, we can read from Eq. ( |37D T>2(R) = — 2e and 
therefore have the third relation that is needed to solve all the three unknown coefficients. 
A calculation leads to 

ai(R) = -2eR/A5 , a 2 = -eR/18, a 3 = . (41) 



The choice of the structure function S^-R) leads to Eq.(p9|) in which the odd parity com- 
ponents disappear, leaving a^(R) undetermined. Another correlation function is needed in 
order to remedy the situation. Since the helicity is u-[Vxu] , we seek a correlation function 
which is related to the helicity of eddies of scale of R: 

T 2 {R) = ([u(R + r,t)-u(r,t)] 

•[Vxu(r + R,t) - Vxu(r,t)]) . (42) 

The proper choice of this correlation function is the crucial idea here. The rest is a straight- 
forward calculation. Using the Navier-Stokes equations to compute the rate of change of 
this quantity we find 

= = -G 2 {R) -2h- vV 2 T 2 {R) , (43) 
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which is the analog of (|37|), and where 

G 2 (R) = {(u(r,t) ■ [V R x [u(r + R,t) x [V R (44) 
x u(r + R, £)]]])} + {term R^-R} . 

The conservation of helicity allows the extraction of two derivatives outside the correlation 



functions. The result can be expressed in terms of our definition (2S): 



G2(R) = ^-^e aX ,e^e VK ,[J^(R) + J^ 7 (-R)] . (45) 



Substituting Eq. (|33|) we find 

G2{R) = 2 ^dR^ HR)[kXkK ~ 5XK] ' (46) 
which is the analog of Eq. (flUI). Using Eq. ( |43D in the inertial interval in the limit v — > 



we find the differential equation 

d 2 a A (R) 1 db 2 (R) 3a 4 (R) h 



dR 2 +5 R dR + R 2 2 ' (47) 

The general solution of this equation is 

a 4 (R) = -hR 2 /30 + «ii?" 3 + a 2 R~ 1 . (48) 

Requiring finite solutions in the limit R — > means that a\ = a 2 = 0. Accordingly we end 
up with Eq. (^). We restate again, that in the preceding derivation, we did not assume that 
our turbulence were isotropic. Other terms of J Q,/37 (R) with j ^ can possibly contribute 
to the total J Q,/37 (R). However, under the assumptions of homogeneity and finite energy 
and helicity dissipations, the R dependence of the isotropic part of J Q,/37 (R) is necessarily 
as stated in (pop. 

It should be noted, that a parallel calculation can be easily carried out for the j > 
sectors of the third order correlation function J Q,/37 (R). In these sectors however, there are 
more irreducible representations than there are in the j = sector (to be exact, there is a 
total of 18 representations for each j > 2. 10 of them with (— ) J+1 parity and 8 of them with 

22 



(— y parity). As a result, for j > we get more unknown functions a q j m (R) than equations, 
and hence we cannot obtain a full solution. Our failure in obtaining a complete set of 
equations for j > sectors, can be attributed to the inadequacy of a turbulence theory that 
involves only moments which are simultaneous in time. The j = sector is indeed unique in 
the sense that its low number of irreducible representations makes a full solution possible. 



V. EXAMPLE: ANALYSIS OF ANISOTROPY IN ATMOSPHERIC 

TURBULENCE 

In this section we present experimental evidence for the utility of the expansion in 
terms of the irreducible representations. The analysis of the experimental data was done in 
collaboration with B. Dhruva, S. Kurien and K.R. Sreenivasan, and the reader is referred 



to the details of 15 . In that work we focused on the 2nd rank tensor structure functions of 



velocity differences 

S a(S (R) = ([u a (r + R) -u a (r)] [t/(r+R) -u^r)}) . (49) 



where homogeneity of the flow is assumed, but not isotropy. This object is symmetric in its 
indices and has even parity in R. In addition it is expected to scale with R in the inertial 
range, with possibly different scaling exponent characterizing contributions of different j. 
We demonstrated in that one can usefully represent S a ^(TV) in the form: 



S a ?(R) = £ a qjm \R\^ B a J m {K) (50) 

qjm 



where a q j m are some numerical coefficients, B^ m (R) are the tensor basis of Sf with a 
definite j, m, and £2 are the exponents associated with the j's irreducible representation. 
The isotropic exponent, (2°\ wm be referred to shortly as £2- We note that the coefficients 



"qjm 



are not arbitrary numerical coefficients, because of the constraints imposed by the 



incompressibility of the flow. In Appendix |A] we derive the explicit form of B q j m (R) and the 
necessary relations among the a<jj m 's. The theoretical development of Appendix |A| serves 
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as a basis for the data analysis; we leave it in the appendix since it is somewhat lengthy. 
Nevertheless the interested reader may find it useful for situations that differ from the one 
treated below. 

The data that we want to consider were taken at Taylor microscale Reynolds numbers 
of about 10,000 JT5| . The data were acquired simultaneously from two single- wire probes 



separated by A = 55 cm nominally orthogonal to the mean wind direction. The two probes 
were mounted at a height of 6 m over a flat desert with a long fetch. The Kolmogorov scale 



was about 0.75 cm. Details of the experimental setup can be found in ref. JL5[]. In that 
reference one can find details of another data set that was analyzed in the same fashion, 
leading to results in agreement with those reviewed here. 

Firstly we tested the isotropy of the flow for separations of the order of A. Using the 
standard Taylor hypothesis, define the "transverse" structure function across A as SV(A) = 
{[ui(Ut) — U2(Ut)} 2 ) and the "longitudinal" structure function as Sl(A) = ([ui(Ut + Ut^) — 
UiiUt)] 2 ) where t& = A/U. If the flow were isotropic we would expect [|T^ 



S T (A) = S L (A) + |^ . (51) 

In the isotropic state both components scale with the same exponent, St,l{A) oc A'' 2 , and 
their ratio is computed from ( |5l]) to be 1 + £2/2 ~ 1-35 where (2 ~ 0.69 (see below). The 
experimental ratio was found to be 1.32, indicating that the anisotropy at the scale A is 
small. We expect that the effects of anisotropy should be most pronounced on the larger 
scales. 

Next, we found the functional form of the basis tensors B^ m (R,) and the algebraic 
relations among the coefficients a q j m according to the discussion in the last paragraph of 
Appendix 

Since the anisotropics are not huge, we focused on the lowest order corrections to the 
isotropic (J = 0) contribution. In other words, we wrote 

S Q/3 (R) = Sg)(R) + S&(R) + Sf 2 (R) + tff 3 (R) . (52) 



24 



We defined the coordinate system such that the mean wind direction was along the 3- 
axis, and the separation between the two probes was along the 1-axis. By assuming axial 



symmetry along the mean wind direction, the tensors S"£q(TL), S'"f 1 (R), S^ 2 (R) were to 
contain only the m = components. In addition, since the two probes measured the velocity 
field only in the mean wind direction, we had only the values of S' 33 (R) in the 1 — 3 plane. 
In such a case, it turns out that only the even j's have a non- vanishing contribution. We 
therefore used the trial tensor: 

fi^(R) = tff (R) + Sfi 2 (R) , 



?a(3 



-taf3 



s°£o(R) 



c » 1 J) 



<2 



(2 + C 2 )5 a/3 -C 2 



Sg,(R) = aSfi 2>q=1 (R) + bS^ q=2 (R). 



Where S^ q=l (R), Sfi 2 , q=2 (R) are given by: 



<S& 9 =i(R) =(i) - 2)5<* - ^ (C^ + 6)5 



A?) 
( »2 



a; 

^(2)12 



(2) 



ro/3 



+2d 2) (cr-2) 



i? 2 



(2) R a RP(k-Rf 



R A 



+ ([Cf] 2 + 3d 2) + 6)^ 



d 2) (d 2) -2) 
i? 2 



(R a kP + R l) k a ){k ■ R) 



(53) 



(2) 



^ , 2(R) = ( [ - (d 2) + 3)(cr ; + 2)^(* • H) a + (cr - 2)-^ 



.(2) 



.(2) 



+ (d 2) + 3)(d 2) + 2)*^ + (2d 2) + l)(d 2) " 2) 



.(2) 



(2) 



.(2) 



X 



R a RP(k ■ R)< 
R~ 4 



ac 



(2)12 



4)(R a kP + R p k a )(k ■ R) 



The vector k was taken to be along the mean wind direction. ( 2 is the isotropic exponent, 

(2) 

while Q is the j — 2 exponent. Co, a, b are the non-universal weights of the components of 
S a ^(R). In order to reduce the number of unknown quantities, the exponent of the isotropic 
part of S a P(R) was assumed to be known: £2 = 0.69. The values of c , a, b, (jp were to be 
found from the experimental data. 

Using spherical coordinates, the trial tensor S* 33 in the 1 — 3 plane took the following 
form: 
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S 33 (R, 9, (f> = 0) = Sf =Q (R, 9,cf> = 0) + S™ 2 (R, 6,(f> = 0) (54) 

2 + 0.69 - 0.69 cos 2 9 



'•<> ( 



a; 



(d 2) + 2) 2 -d 2) (3d 2) + 2)cos^ 



+2d 2) (d 2) -2)cos 4 ^ 

\" J r,*(2) , oV/1 2) , :m , i2i,,, 1.2, ,._2, 



+b[~) [(cr ; + 2) (cr ' + 3) - cr (scr j + 4) ^ e 



+(2cr + i)(cr-2)cos^ 

where 9 is the angle between k and R. The fitting of the trial tensor to the data was done 
along two paths in the (R, 9) space: 

• 9 = A single probe measurement. 

• Rsin9 = A Two probes measurement. 

Fig. 1 shows the best fit to the data. For each type of data, two fits were performed: 
A fit of the isotropic part only (panel a), and a fit of both isotropic and j = 2 components 
(panel b). The excellent fits in panel (b) is a good support for the present mode of analysis. 

In Ref. Jl5| Fig. 3 we showed the determination of from a least-square fit. The optimal 



(2) 

value of this exponent and the uncertainty determined from that plot is Q — 1.38 ± 0.15. 
It should be understood that the exponent Q (and also Q that is unavailable from the 
present measurements) are just the smallest exponents in the hierarchy that characterizes 
higher order irreducible representations indexed by j. The study of these exponents has only 
begun here, and considerable experimental and theoretical effort is needed to reach firm 
conclusions regarding their universality and numerical values. We expect the exponents 
to be a non- decreasing function of j, explaining why the highest values of j are being 
peeled off quickly when R decreases. Nevertheless, the lower order values of can be 
measured and computed. In Ref. 113 we presented an additional set of experimental data, 



(2) 

and demonstrated that the numerical value of Q appears universal. 
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VI. CONCLUDING REMARKS 



The aim of this paper has been to introduce the exploration of the scaling properties 
of turbulent statistics in the anisotropic sectors. The main novel theoretical development 
is described in Sect. 2. It is explained there that the linearity of the equations for the fully 
unfused correlation functions together with the invariance to rotations, foliates the solutions 
into sectors characterized by the j, m designations of the irreducible representations of the 
SO(3) symmetry group. As a consequence we expect the different sectors to be characterized 
by different (anomalous) scaling exponents. This observation opens up a new and interesting 
research arena for theory and experiments. 

In Sect. 3 we presented a derivation of the form of the second order structure function 
in the higher j-sectors, and used a recent experimental data analysis as an example of the 
utility and importance of the present approach. The main result of this section, besides 
the theoretical forms that can be used for further data analysis, is that the scaling range in 
turbulence is much larger than expected. One just needs to acknowledge the existence of 
anisotropic contributions to obtain scaling ranges that go all the way from the Kolmogorov 
scale to the outer scale of turbulence. 

It is our belief that this paper does not exhaust the issue of anisotropic contributions 
to turbulent statistics. It is only the beginning of a rich research program that should be 
carried simultaneously by experimentalist and theorists. 
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APPENDIX A: THE GENERAL FORM OF THE 2ND RANK TENSOR 

In this appendix we discuss the general structure of the second rank correlation functions 

F a/3 (R) = (u a (r + R)v?(t)) , (Al) 

In (|A1|) homogeneity of the flow is assumed, but not isotropy. Note that this object is more 
general than the structure function S a/3 in being nonsymmetric in the indices, and having no 
definite parity. In light of the discussion in Sect. 2, when we expand this objects in terms of 
tensors with definite j, m, we expect each component to have a distinct dependence on the 
amplitude R = |R|. Accordingly, we wish to find the basis functions B q ^ m (R), with which 
we can represent F a/3 (R) in the form: 

F^(R) = £ a qjm (R)B^ m (R) (A2) 

qjm 

and derive some constraints among the functions a q j m (R) that result from incompressibility. 
We shall see, that due to the isotropy of the incompressibility conditions, the constraints 
are among a q j m (R) with the same j,m only. 

We begin by analyzing the incompressibility condition: An incompressible flow with 
constant density is characterized by the relation: 

d a u a (r,t) = 

as a result, one immediately gets the following constraints on F al3 (R): 

d a F a/3 {R) = 
d p F a0 (R) = 0. 

Plugging the trial tensor ( |A^ ) into the last two equations we obtain 2 equations connecting 
the different a Q ,- m : 
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d a £ a qjm (R)Bgf m (R) = (A3) 

qjm 

a^a cm (i?)5j(R) = 

qjm 

We first notice that the differentiation action is isotropic. As a result, if T Q/3 (R) is some 
arbitrary tensor with a definite j, m transformation properties, then the tensor <9 Q T Q/3 (R) will 
have the same j, m transformation properties. Components with different j, m are linearly 
independent. Therefore equations flA3| ) should hold for each j, m separately. 

Next, we observe that ( |A3|) are invariant under the transformation F a @ — > F 13 ". As a 
result, the symmetric and anti-symmetric parts of F af3 should satisfy (|A3|) independently. 
To see that, let us write F al3 as a sum of a symmetric term and an anti-symmetric term: 

F aP = + F a0 we then get . 

daF* = d a Ff + d a Ff = d a F§ a - d a F? a = 
dpF * = dpFf + dpFf = 

from which we conclude: 

d a Ff = d a Ff = 

Finally, (|A3|) is invariant under the transformation F a/3 (R) — ► F a ^{— R) and as a 
result the odd parity and the even parity parts of F a/3 should fulfill (|A3[) independently. We 
conclude that a necessary and sufficient condition for ( |A3|) to hold is that it holds separately 
for parts with definite j, m, definite symmetry in the a, (3 indices and a definite parity in R: 

„ „ summation is over 

5«E a ^m(|R|)S^ TO (R) = 

i Bqj m with definite symmetries 

where the summation is over q such that B^ m has a definite indices symmetry and a definite 
parity. 

According to fllPl) we can write these B^ m as: 
1. (— Y parity, symmetric tensors: 

29 



. B"j m (R) = R~ j 5 a/3 $ jm (R), 

• Sg m (R) = iT' [i^ + i^<9 a ] $, m (i? ), 
. sg m (R) = i2-J- 2 if*i^$ im (i2), 

. B^ m {R)^R-^d^ jm {R)- 

2. (— parity, anti-symmetric tensors: 

3. (— y +1 parity, symmetric tensors: 

• Sg m (R) = R^- 1 [R a e^R,d u + i^^ct] $, m (i?), 

• Bf jm (R) = R^ +1 [e^R^d a + e^R^d?] $ jm (R). 

4. (— Y +1 parity, anti-symmetric tensors: 

. B2? m (R) = R-i-^R^UR), 
. Bf jm {R) = R-^e^d^ jm {R). 

In order to differentiate these expressions we can use the following identities: 

R a d a R% m (R) = (R c Y jm (R) , 

d a d a R% m {K) = [C (C + 1) - j(j + 1)] R^Yjmix) 



which give rise to: 



R a d a <S> jm {R)=j<S> jm {K) . 

d a d a ^ ]m {K) = o . 

From this point, it is a matter of simple (though somewhat lengthy) algebra to derive the 
differential constraints among a q j m (R). The results are as follows: 
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1. q G {1,7,9,5} 

a'ijm(R) ~ jR'^ijm + ja' 7Jm - 3 2 R- X a^ jm + a' 9 jm + 2R~ l a^ m = , (A4) 
R~ l a\,jm + a' 7 jm + 3i? _1 a 7Jm + (j - 1) a' b jm - (j 2 - 3j + 2) i? _1 a 5Jm = . 

2. q e {3} 

a'z,im - jR^asjm = , (A5) 

a 3,jm + R la 3,jm = . 

Notice: These equations have no solutions other than: a 3jm (i?) = 0. 

3. qe{8,6} 

*8, 3m + 3irV, jm + (j - l)a' 6Jm - (f - 2j + l) R- l a,, jm = . (A6) 

4. gG{4,2} 

R^a^jm - a' 2 jm + (j - l)R~ 1 a 2 , jm = . (A7) 

There are obviously more unknowns than equations, since we merely exploited the incom- 
pressibility conditions. Nevertheless, we believe that the missing equations that arise from 
the dynamical hierarchy of equations will preserve the distinction between a g j m of different 
j,m (again, due to the isotropy of these equations). 

Note also, that the above analysis holds also for the second-order structure function 

S afS (R) = ([u a {r + R) - u a (r)] f^(r + R) - u p {r) 



Only that in this case we should only consider the representations q — 1, 7, 9, 5 for even j and 
the representations q = 8, 6 for odd j. This follows from the fact that S^^R) is symmetric 
with respect to its indices and it has an even parity in R. Also, in that case, it is possible 
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to go one step further by assuming a specific functional form for the a q j m (R). We know 
that the S a/3 (R,) is expected scale in the inertial range, and we therefore may assume: 

where c q j m are just numerical constants. If we now substitute this definition into the equa- 
tions QA41JA6D, we get a set of linear equations among the c q ^ m . These relations can be easily 



solved and give us two possible tensors for even j (q = 1, 7, 9, 5) and one tensor form for odd 
j (from q = 8,6). This kind of approach was taken in the two-probes experiment which is 
described in Sect. 0. 
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FIG. 1. The structure functions S 33 for = and for non-zero computed for set I. The dots 

are for experimental data and the line is the analytic fit. Panel (a) presents fits to the j = 

component only, and panel (b) to components j = and j = 2 together. 
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